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ABSTRACT 


The  performance  potential  of  a  room 
temperature  magnetic  heat  pump  utilizing 
Gadolinium  and  operating  on  an  Ericsson  Cycle  was 
investigated  at  magnetic  flux  densities  of  2  and 
7-Tesla  which  represent  the  upper  limits  of 
conventional  and  superconducting 
electromagnetics,  respectively.  At  a  coefficient 
of  performance  of  5,  a  7 -Tesla  system  would 
provide  a  cooling  capacity  of  at  best  1200  BTU 
per  hour  per  pound  of  Gadolinium  while  a  2-Tesla 
system  would  operate  at  approximately  130  BTU  per 
hour  per  pound  of  Gadolinium.  Magnetic  circuit 
efficiency  was  not  determined  but  must  be  high 
(95  -  percent  or  better)  in  order  for  the 
magnetic  heat  pump  performance  to  compete  with 
conventional  cooling  systems.  It  is  unlikely  the 
magnetic  heat  pump  investigated  could  approach 
the  performance  and  compactness  of  the 
conventional  cooling  systems  unless  field 
strengths  much  greater  than  7-Tesla  are  possible. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  Center's  IR/IED  project  entitled 
"Magneto/ Thermal  Heat  Pump",  Work  Unit  2723-135,  Program  Element  62766N,  Task  Area 
ZF61512001.  This  work  was  performed  in  the  Engines  Branch,  Power  System  Division  of 
the  Propulsion  and  Auxiliary  Systems  Department. 


INTRODUCTION 

Magnetic  heat  pumps  are  based  on  the  magnetocaloric  effect  by  which  certain 
materials  exhibit  an  increase  (or  decrease)  in  temperature  with  the  application  (or 
removal)  of  a  magnetic  field.  Magnetic  heat  pumps  originally  found  applications  only 
at  temperature  near  absolute  zero  using  paramagnetic  materials  Using  a 

ferromagnetic  material  having  a  Curie  point  in  the  ambient  temperature  range,  a 
magnetic  heat  pump  for  air  conditioning  is  possible  by  utilizing  a  suitable 
regenerative  thermodynamic  cycle  ^)»  (3)t  Rare  earth  materials  are  effective 
in  this  application  and  Gadolinium  (GD)  with  a  Curie  point  of  68  °F  (293  °K) 
is  a  good  working  material.  At  Its  Curie  point,  the  application  of  a  7-Tesla 
field  to  GD  causes  a  temperature  rise  of  25°F  (14  °K)  or  the  release  of  1.7 
BTU/lb  (4  KJ/Kg)  of  heat  under  isothermal  conditions. 
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Magnetic  heat  pumps  are  of  interest  for  several  reasons.  The  number  of 
moving  parts  in  the  system  should  be  small  and  thus  high  reliability  is 
expected.  For  a  ferromagnetic  material  like  Gadolinium,  an  Ericsson  or 
Stirling  cycle  are  similar  and  cycle  efficiency  can  approach  the  Carnot  cycle 
efficiency.  This  should  provide  considerable  improvement  over  conventional 
room  temperature  refrigeration  systems  which  generally  operate  on  a  Rankine 
type  cycle. 

In  this  study,  a  magnetic  heat  pimp  operating  on  an  Ericsson  cycle  was 
evaluated  parametrically  to  determine  expected  thermodynamic  performance. 
Losses  associated  with  the  generation  of  the  proper  magnetic  field  are  not 
determined  in  this  report  however  the  effect  of  the  magnet  circuit  efficiency 
on  the  cycle  performance  is  shown.  A  2  and  7  -  Tesla  field  strength  were 
examined  with  2-Tesla  representing  the  limit  of  an  electromagnet  and  7-Tesla 
the  limit  of  a  superconducting  magnet. 

BACKGROUND 

A  typical  magnetic  heat  pump  consists  of  a  porous  mass  of  ferromagnetic 
material,  Gadolinium  plates  in  this  case,  located  in  a  vertical  column  of 
fluid  with  an  energizing  magnet  surrounding  the  column  at  the  location  of  the 
Gadolinium  plates.  At  least  one  of  the  working  substances  must  move.  One 
case  is  to  allow  the  column  of  fluid  to  remain  stationary  and  to  move  the 
Gadolinium  located  in  the  column  up  and  down.  This  also  requires  moving  the 
magnet  which  surrounds  the  column  with  the  Gadolinium,  resembling  test  setups 
by  Brown  ^  .  One  of  the  dravAacks  of  this  setup  is  that  the  magnet  and 

Gadolinium  are  separated  by  the  column  wall  plus  gaps  to  allow  for  movement, 
all  resulting  in  a  degradition  of  magnet  performance. 

The  arrangement  investigated  keeps  the  Gadolinium  (and  thus  the  magnet) 
stationary  and  moves  the  column  of  fluid  back  and  forth  through  the 
Gadolinium.  This  allows  the  Gadolinium  to  be  in  as  close  a  proximity  to  the 
magnet  as  possible  to  minimize  magnetic  losses.  Moving  the  liquid  column 
requires  the  system  to  be  about  twice  as  long  as  the  above  stationary  fluid 
case  using  the  same  Gadolinium  geometry. 

Figure  1  shows  a  schematic  of  the  magnetic  Ericsson  cycle  investigated. 
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Figure  2  and  3  show  the  magnetic  Ericsson  cycle  entropy- temperature  diagram 
and  the  fluid  column  and  Gadolinium  temperature  profiles,  respectively.  The 
working  fluid  was  assumed  to  be  water  and  all  points  In  the  cycle  are  numbered 
similarly  In  the  figures.  The  column  would  be  mounted  vertically  with  the 
cold  end  at  the  bottom  and  the  hot  end  at  the  top  to  avoid  fluid  mixing  caused 
by  density  differences  in  the  fluid  column.  Heat  exchangers  are  located  at 
the  top  of  the  column  to  remove  the  heat  of  magnetization  from  the  fluid  and 
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Figure  1:  Magnetic  Ericsson  Cycle  Heat  Pump  Schematic 
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Figure  2:  Magnetic  Ericsson  Cycle  Entropy-Temperature  Diagram 
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at  the  bottom  to  interface  wi th  the  refrigeration  load.  Centrally  located  in 
the  column  is  the  stationary  Gadolinium  surrounded  by  its  magnet.  The  fluid 
column  is  composed  of  essentially  three  sections.  The  cold  column,  the 
regenerative  column,  and  the  hot  column  as  shown  in  figure  1  and  3. 

The  Ericsson  cycle  ideally  consists  of  two  isothermal  processes  and  two 
constant  pressure  processes.  Referring  to  figure  1  to  3,  the  cycle  operation 
is  as  described  below: 

1- 2:  Ideally,  this  is  an  Isothermal  process  during  which  the  magnetic 

field  is  decreased  from  some  maximum  value  to  some  minimum  value. 

The  demagnization  produce  a  cooling  effect  in  the  Gadolinium  plates 
and  heat  is  transferred  from  the  water  flowing  thru  the  Gadolinium 
plates  to  maintain  the  plates  at  constant  temperature  while  chilling 
the  water.  Because  temperature  differences  are  required  between  the 
water  and  the  plates  for  heat  transfer  to  actually  occur, 
isothermal  operation  is  possible  from  state  points  1'  to  2  only  with 
1  to  1'  representing  a  loss  due  to  the  finite  temperature 
differences  required  for  heat  transfer. 

2- 3:  This  Is  a  regenerative  process  during  which  heat  is  transferred  from 

the  water  column  to  heat  the  plates.  Because  entropy  change  in  the 
Gadolinium  at  low  or  no  field  shows  a  definite  peak  centered  on  the 
Curie  temperature,  an  appropriate  variation  of  the  applied  magnetic 
field  during  this  regeneration  process  is  necessary  so  that  the 
entropy  change  from  statepoints  2  to  3  can  be  balanced  with  an  equal 
but  opposite  entropy  change  during  regenerative  process  5-6. 

3- 4:  The  hot  water  column  passes  through  the  Gadolinium  plates  and  acts 

as  a  continuation  of  the  regeneration  process  2-3.  Heat  is 
transferred  from  the  fluid  to  the  plates. 

4- 5:  Ideally  this  Is  an  isothermal  process  during  which  the  magnetic 

field  Is  Increased  to  its  maximum  value.  The  magnization  of  the 
Gadolinium  produces  a  heating  effect  and  this  heat  is  transferred  to 
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the  water  flowing  through  the  plates  to  maintain  the  plate  at 
constant  temperature.  As  in  process  1-2,  temperature  differences 
are  required  between  the  water  and  the  plates  for  heat  transfer  to 
actually  occur.  Thus  isothermal  operation  is  only  possible  from 
state  points  4'  to  5  with  4  to  4*  representing  a  loss  due  to  the 
establishment  of  the  finite  temperature  differences. 

5- 6:  This  is  a  regenerative  process  during  which  heat  is  transferred  from 

the  Gadolinium  plates  to  the  water  flowing  through  them.  The 
magnetic  field  is  at  its  maximum  and  held  constant  while  the  plates 
are  cooled  in  the  regeration. 

6- 1:  The  cold  column  passes  through  the  Gadolinium  plates  and  acts  as  a 

continuation  of  the  regeneration  of  process  5-6.  Heat  is 
transferred  from  the  plates  to  the  fluid. 

The  fluid  is  pumped  toward  the  cold  end  heat  exchanger  for  cycle  points 
1-2 -3-4.  Then  the  direction  of  pumping  is  reversed  toward  the  hot  end  heat 
exchanger  for  cycle  points  4-5-6-1.  When  the  cold  column  reaches  the  cold  end 
heat  exchanger  or  the  hot  column  reaches  the  hot  end  exchanger,  external  loops 
exchange  the  cooling  load  and  the  heat  of  magnetization,  respectively. 


CYCLE  DETAILS 


To  simplify  the  analysis  it  was  assumed  that  the  Gadolinium  geometry 
would  be  parallel  flat  plates  as  shown  in  figure  4.  Water  flowing  between  the 
plates  either  adds  heat  to  or  removes  heat  from  the  plates.  Since  the  plates 
act  as  a  restriction  to  the  flow  in  the  assumed  constant  cross-sectional  area 


fluid  column  channel,  the  water  velocity  between  the  plates,  Vp,  is  greater 


than  the  fluid  column  velocity,  V^,  and  is  given  by. 


(1) 
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MAGNETIC 

FIELD 


~~  """"Sp 


WCL  -  COLUMN  WIDTH 

Lp  -  PLATE  LENGTH 

Hch  =  CHANNEL  HEIGHT 

tp  -  PLATE  THICKNESS 

Sp  -  PLATE  SPACING 

Vp  *  PLATE  CHANNEL  VELOCITY 


Figure  4:  Gadolinium  Plate  Geometry 


DTNSRDC  PAS-82-11 


Velocity  through  the  plates  Is  Important  since  it  helps  determine  whether  the 
water  side  heat  transfer  coefficient  is  based  on  a  laminar  or  turbulent  heat 
transfer  correlation.  Laminar  flow  heat- transfer  coefficients  were  based  on 
average  conditions  from  tabulated  data  in  reference  4  for  flat  ducts. 
Turbulent  heat  transfer  coefficients  were  obtained  using  the  following 
correlation  for  turbulent  flow^: 


hD  =  0.023  Re  °*8  Pr  0,33 
H 


(2) 


The  overall  transmittance  between  the  Gadolinium  plates  and  the  water  was 
calculated  by 


U  = 


(3) 


Table  1  presents  physical  properties  of  Gadolinium  that  were  used  in  the 
calculations  in  this  report.  For  all  calculations  in  this  report  average 
conditions  are  used.  Thermodynamic  properties  of  Gadolinium  in  the  presence 
of  a  magnetic  field  were  obtained  from  reference  6  and  are  plotted  in  figure  5 
as  a  function  of  entropy  versus  temperature  for  magnetic  fields  from  0  to  7 
Tesla.  This  data  was  tabulated  for  use  in  the  computer  program  (See  Appendix 
A)  that  was  written  to  perform  the  cycle  calculations.  Figure  5  shows  the 
nonlinearity  of  the  entropy  change  at  low  field  strength  around  the  Curie 
temperature  (  293  °K).  This  nonlinearity  affects  the  cycle  regeneration 
processes  and  necessitates  the  variation  of  the  field  strength  discussed 
earlier  in  order  to  balance  the  entropy  change  during  the  regeneration  process 
between  the  low  field  strength  state  and  high  field  strength  state. 
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Table  1  -  Physical  Properties  of  Gadollnuim  at  68°F  (20°C) 


Molecular  Weight  (Mqq):  157.25  Ibm/mole  (346.74  Kg/mole) 

Oenslty  (  pgq):  490.7  lbm/ft3  (7.86  gm/cm3) 

Spelclfc  Heat  (cvQD)  0.071  BTU/lbm  °F  (.71  CAL/gm  °C) 

Thermal  Conductivity  (kQD):  5.8  BTU/Hr  ft  °F  (.1  Watt/cm  °C) 

In  evaluating  conditions  around  the  cycle,  one  of  the  losses  that  affects 
performance  Is  frictional  flow  losses  between  the  plates.  It  was  assumed  that 
frictional  losses  were  dissipated  as  heat  which  heated  the  working  fluid.  The 
amount  of  heat  generated  by  friction  Is  found  by 


QF  =  4.6263  H]  y  Vp.  t  .  N  p 

(4) 

where 

L 

H1  =  Head  Loss  *  ~T~  *  IT-  *  g 

H 

(5a) 

and 

A«.  *  Plate  Surface  Area  *  2  .  Lp  .  HCH 

(5b) 

In  the  laminar  flow  region,  friction  factor  can  be  found  by  assuming  the  plate 
geometries  satisfy  conditions  for  two-dimensional,  steady  flow  between 
infinite  parallel  plates.  Then  solving  the  Navler-Stokes  equation  to  obtain 
heat  loss  and  equating  to  equation  5  yields 


96 


Re  *  (1  +  J-  ) 
hCH 


(6) 


In  the  turbulent  flow  region,  friction  factor  between  the  plates  was  asstmed 
to  be  comparable  to  the  friction  factor  for  a  smooth  pipe  In  turbulent  How 
and  was  estimated  from  the  following  curve  fit  of  the  smooth  pipe  friction 
factor  curve, 
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(7) 


No  friction  factor  was  determined  in  the  transitional  flow  region  and  the  flow 
was  assumed  to  go  from  laminar  to  turbulent  at  a  Reynolds  number  (Re)  of  3000. 

Mixing  losses  can  occur  in  the  fluid  column  and  at  the  plate  inlets  and 
exits.  Ideally,  a  horizontal  cut  through  the  fluid  column  would  show  the 
entire  layer  of  fluid  to  be  at  uniform  temperature  and  velocity.  However, 
moving  the  fluid  column  up  and  down  can  result  in  velocity  profile  being 
established  across  the  column  due  to  the  presence  of  the  stationary  walls. 
Developing  turbulent  flow  has  a  nonuniform  velocity  profile  which  essentially 
becomes  uniform  when  fully  developed,  except  close  to  the  walls.  Laminar  flow 
has  a  parabolic  velocity  profile  when  fully  developed.  The  velocity  profile 
causes  fluid  of  one  temperature  to  be  introduced  into  a  region  of  fluid  at 
another  temperature  causing  mixing  losses  due  to  the  temperature 
differential.  In  the  analysis,  mixing  losses  were  considered  to  cause  a 
temperature  differential  across  the  hot  and  cold  column  (  ATmi.x). 

Other  loses  can  occur  in  the  system  but  were  not  specifically 
considered.  Other  losses  can  include: 


o  Inlet  and  exit  losses  at  the  Gadolinium  plates 
o  Friction  between  the  water  column  and  the  container  walls 
o  Pumping  power  to  overcome  losses 

Although  these  losses  are  not  directly  considered,  the  first  two  could  be 
considered  to  be  grouped  in  with  the  mixing  loss.  Pumping  power  will  reduce 
the  system  overall  performance  but  would  not  change  the  cycle  state  points 
provided  the  losses  requiring  pumping  power  are  considered.  Any  losses  over 
and  above  those  considered  in  this  analysis  would  further  reduce  system 
overall  performance.  The  loses  considered  in  this  analysis  should  give  an 
overall  indication  of  the  system  performance. 

CYCLE  ANALYSIS 

In  the  analysis,  a  cold  end,  external  sink  temperature,  Tcolc|,  that  is  to 
be  maintained  by  the  heat  pump  is  specified,  as  is  the  hot  end,  external  sink 
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temperature,  Thot.  In  order  for  heat  transfer  to  occur  between  the  heat  pump 
fluid  and  the  external  heat  exchanger  fluids,  some  heat  exchanger  temperature 
difference,  TH£,  is  required  and  is  specified  in  this  analysis  and  assumed  to 
be  the  same  for  both  the  hot  and  cold  end  heat  exchangers.  Thus  cold  and  hot 
columns  must  be  operated  at  some  lower  and  higher  temperatures,  respectively, 
tu  effect  the  heat  transfer  to  the  external  heat  exchangers  and  are  define  by, 


T(c)  3  Tcold  "  aTHE 

(8) 

T(H)  3  Thot  +  aTHE 

(9) 

T(C)  is  the  temperature  of  the  cold  fluid  column  between  state  points  f  and  1 
and  T(H)  is  the  hot  fluid  column  temperature  between  state  points  c  and  4  for 
the  ideal  cycle  of  figures  2  and  3.  However,  due  to  friction,  mixing,  and 
heat  transfer  temperature  differences,  process  f  1  and  process  c  4  are  not 
isothermal  and  T(C)  and  T(H)  were  specified  to  be  the  average  cold  and  hot 
column  temperatures  between  state  points  f  and  1  &  c  and  4,  respectively. 

Figure  6  shows  the  actual  cycle  considering  losses  and  the  heat  transfer 
tempeature  differences.  Figure  7  shows  the  actual  fluid  column  and  Gadolinium 
temperature  profiles. 

The  system  was  evaluated  on  a  per  pound  of  Gadolinium  basis.  Cooling 
load  was  evaluated  by, 

«i-2  *  <^4-^  *  459-7)-  -'-rj—  •  <60 “  rpm'  -  V2  (10) 

Where  S2  and  Sj  are  the  total  entropy  of  Gadolinium  at  and  Tj 

respectively.  During  this  process  ( 1 — ►  2 )  the  magnet  field  is  decreased  from 
a  maximum  value  at  T^  to  near  zero  at  Tg  and  values  of  entropy  are  obtained 

from  figure  5.  Since  this  is  a  cooling  process,  fricton  generated  by  the 

water  flowing  between  the  plates  reduces  the  cooling  load  by  which  is 

evaluated  using  equation  4.  Not  all  values  for  equation  10  and  subsequent 
equations  are  known  at  the  start  but  initial  assumptions  were  made  and 
subsequent  iterations  in  the  computer  program  to  evaluate  the  cycle  provided 
rapid  convergence  to  the  correct  values. 
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TEMPERATURE 


Figure  7:  Actual  Column  and  Gadollniun  Temperature  Profiles 
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The  entropy  change  at  the  hot  and  cold  end  will  be  equal  in  magnitude  but 
opposite  in  sign  for  the  ideal  isothermal  case.  Thus  for  the  Gadolinium, 


( S2I  '  S1I}  T(C)  =  (S4I 


S5I)  T(H) 


However  due  to  losses  and  temperature  difference,  the  entropy  change  at  the 
hot  and  cold  end  are  no  longer  necessarily  equal  in  magnitude.  It  was  assumed 
that  the  temperature  -  entropy  line  from  state  points  2  to  4  was  parallel  to 
the  maximum,  constant  field  temperature  -  entropy  line  from  state  points  5  to 
1  in  order  to  balance  the  loads  on  the  regenerative  portion  of  the  cycle. 

This  requires  some  magnetic  field  variation  from  state  points  2  to  4  and 
assuming  the  lines  are  linear,  entropy  at  point  4  for  Gadolinium  is  found  by 


S4  =  S4I  +  (S4I  '  S2I) 


T4  _  T(H) 


Thus  the  heat  load  to  the  hot  column  is  found  from. 


rT4  ♦  T5 


+  459.7)  . 


(S4  -S5) 


.(6°  NRpM  )+  Qp 


with  S5  being  at  the  maximum  field  strength  and  Qp  4_5  being  heat  added  due  to 
friction  between  the  plates,  which  now  adds  to  the  hot  end  heat  rejection 
load.  The  regenerative  heat  loads  are  found  from 


Q2-3  =  60  cvGQ  MRPM  (T3"  V 


Q5-6  =  60  cvGD  NRPM  (W 


Likewise  the  cold  and  hot  columns  pass  through  the  Gadolinium  prior  to  the 
change  in  magnetic  field  and  exchange  heat  with  each  other  where, 
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Q6-l  =  60  CvQD  NRPM  (T6  '  V 


Q3-4  =  60  cv™  NRPM  (T4  "  V 


Heat  transfer  to  and  from  the  water  is  related  to  the  mass  flow  rate  of 
water  and  the  channel  geometry.  Based  on  a  one  pound  weight  of  Gadolinium, 
the  number  of  plates,  Np,  in  the  system  is, 


WT^gd 


Thus  the  column  width,  WCH,  is 


WCH  =  (tP  +  V  NP 


The  free  flow  area  between  the  plates,  App,  is 


Vf  = 


HCHSpNP 


Thus  the  mass  flow  rate  of  water,  mu,  through  the  plates  is, 


rlly  =  3600  P^App  Vp 


In  calculating  temperature  changes  in  the  water,  an  estimate  of  the  time 
for  each  of  the  three  column  sections  {cold,  regenerative,  and  hot)  to  pass 
through  the  plates  was  required,  the  column  passes  through  the  plates  twice 
to  complete  a  cycle.  The  overall  tansmittance  and  temperature  differential  in 
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the  regenerative  column  and  the  hot  and  cold  columns  during  the  change  in 
magnetic  field,  were  close  and  thus  the  ratio  of  the  heat  transfers  in  these 
areas,  provided  a  close  estimate  of  time  for  each  column  section  to  pass 
through  the  plate. 


*T0T 


=  q2-3  +  Q5-6  +  2*  (Ql-2  +  ^  +  q4-5  ’  <V 


4-5 


)  (22) 


then  for  the  cold  column 


cold 


Ql-2  +  Fl-2 

’’tot 


60  N 


RPM 


(23) 


for  the  regenerative  column 

*  Q2-3  +  Q5-6 

REG  ^0~ 

and  for  the  hot  column 


60  N 


RPM 


(24) 


v4-5  - 


4-5 


hot 


60  N 


nOT 


RPM 


(25) 


Determining  the  temperatures  at  the  state  points  is  accomplished  based  on 
a  few  assumptions.  It  is  assumed  at  state  points  1  and  4,  that  the 
temperature  of  the  fluid  and  Gadolinium  are  the  same  prior  to  the  changing  of 
the  magnetic  field  (see  figs  6  &  7).  Also  the  mean  temperature  of  the  cold 
and  hot  fluid  columns  prior  to  changing  the  magnetic  field  was  equal  to  the 
values  defined  by  equations  (8)  and  (9)  respectively  (see  figs.  5  &  7).  Thus, 
for  the  cold  column 


T(C)  = 


(26a) 


and  for  the  hot  column 
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T(H)  =  JL-^—! -  (26b) 

If  no  losses  or  temperature  difference  were  present,  lines  f-1  and  c-4  would 
be  isothermal  as  previously  discussed.  However  the  two  temperatures  are  not 
equal  but  are  different  by 


Tf  -  T1  ■  lTM1x  +  lTF 


5-6 


(27a) 


for  the  cold  end  and  for  the  hot  end 


T4  -  Tc  ’ 


AT, 


Mix 


AT. 


2-3 


(27b) 


The  mixing  temperature  change,  Tffl1x,  was  difficult  to  evaluate  and  was 
assumed  to  be  some  fixed  temperature  difference  across  the  hot  and  cold 
column.  The  fluid  flowing  through  the  plates  generates  friction  which  in 
effect  adds  heat  to  the  fluid.  The  temperature  increase  of  the  fluid  caused 
by  friction  is  obtained  in  the  different  column  sections  by  use  of  ether 
laminar  or  turbulent  equations  (4  to  7)  and 

Qc 

6-1  (28a) 


ATr 


6-1 


*W  V  *001  d 

rw 


ATr 


5-6 


5-6  CPW  t  REG 


(28b) 


AT, 


4-5 


4"5  CP.„  W 


(28c) 


AT. 


3-4 


3-4 


*W  Cpw  W 


(28d) 


AT, 


2-3 


2-3 


cpm  *Reg 


(28e) 
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(28f ) 


AT  *  - ±-± - 

1-2  *W  Cp  *REG 
rw 

In  order  for  the  Gadoliniun  and  the  fluid  to  achieve  the  same  temperature  at 
points  1  and  4  heat  must  be  transfered  between  the  two.  Thus  the  heat  gain  or 
loss  of  the  Gadolinium  was  reflectd  in  the  fluid  as  a  heat  toss  or  gain 
respectively.  The  change  in  temperature  of  the  fluid  in  exchanging  heat  with 
the  Gadolinium  was  at  the  cold  end, 


AT  „ 


-W-GD 


60  NRPM  vGD  (T6  "  V 

\  co  ^old 


and  at  the  hot  end 


(29a) 


AT„ 


w-GD 


60  H  RPM  VGD  (T4~  V 
lhW  Cpw  W 

w 


(29b) 


The  temperature  changes  are  shown  in  figure  8  which  will  be  discussed  shortly. 

Additional  information  is  still  needed  to  obtain  Tg  and  T^  for  equation 
29.  This  involves  finding  the  average  temperature  difference  in  the 
regenerator  between  the  plates  and  the  fluid.  Thus, 


AT 


*5-6 


REG5-6  U5-6  AS  *REG 


and 


AT, 


REG 


2-3 


t__2-3__ 

J2-3  AS  *REG 


(30a) 


(30b) 
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Figure  8:  Temperature  Changes  in  Cold  and  Hot  Column 


with  the  assumption  that, 
T6  -  Tf  =  aTREG 


Ola) 


Y  T3  =  aTREG 


(31b) 


The  average  temperature  change  of  the  cold  and  hot  fluid  column  due  to 
the  change  in  magnetic  field  is  found  by 


\  CPW  Yld 


(32a) 


and  for  the  hot  end 


H  ^  CPW  W 


(32b) 


The  temperature  difference  between  the  Gadolinium  and  the  water  to  effect  the 
heat  transfer  in  the  cold  and  hot  column  during  field  change  is  found  by 


0 

0-2  +  0-2 


Cl'-2  Ul-2  AS  Yld 


(33a) 


and  for  the  hot  column 


r  -°4-5  ‘Vs 

H4'-5  U4-5AS  hot 


(33b) 
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All  the  above  temperatures  are  Illustrated  In  Figure  6. 


All  the  heat  loads  and  temperatures  have  been  described  up  to  this 
point.  Repeated  iteration  with  the  equations  presented  led  to  a  rapid 
solution  of  the  cycle  in  the  computer  program  that  was  written  to  evaluate  the 
cycle  (see  Appendix  A).  All  that  remains  is  to  predict  performance.  The 
cooling  load  that  can  be  transferred  at  the  cold  end  to  the  external  load  is 
less  than  the  change  of  magnet  energy  in  the  Gadollniun  due  to  friction  and 
heat  transfer.  Likewise,  friction  and  heat  transfer  between  the  Gadolinium 
and  fluid  make  the  heat  rejected  at  the  hot  end  differ  from  the  magnetic 
energy  change.  Figure  8  shows  in  more  detail  the  variation  of  Gadolinium  and 
fluid  temperature  in  the  cold  and  hot  column.  Line  If  represents  the 
temperature  distribution  in  the  cold  column  when  the  magnetic  field  starts  to 
decrease  at  point  1.  Line  ab  is  the  temperature  distribution  in  the  cold 
column  fluid  after  the  drop  in  magnetic  field  is  completed  and  the  fluid  has 
been  chilled  by  transferring  heat  to  the  Gadolinium.  Thus  area  labf 
represents  the  amount  that  the  fluid  is  chilled  and  is  equal  to  the  value 
obtained  by  equation  (10).  However  only  the  area  a'abf'  can  be  used  to  chill 
the  external  loop.  Part  of  the  cooling  load,  area  la'f’f  must  be  kept  in  the 
loop  to  offset  frictional  heating  of  the  fluid  and  heat  transfer  from  the 
Gadolinium  to  the  fluid  during  process  6-1.  Line  a’f’  represents  the 
temperature  profile  in  the  cold  column  after  heat  transfer  is  completed  with 
the  cold  end  external  loop. 

In  the  hot  end,  line  c4  is  the  temperature  distribution  in  the  hot  column 
at  the  start  of  the  increase  in  magnetic  field  at  point  4.  Line  de  is  the 
temperature  distribution  in  the  fluid  after  the  magnetic  field  is  brought  to 
Its  maximum  value  and  heat  generated  in  the  Gadolinium  has  been  transferred  to 
the  fluid.  Area  ec4d  is  the  heat  generated  In  the  Gadoliniun  from  equation 
(13)  that  is  transferred  to  the  fluid.  However  area  ec'd'd  is  the  heat  that 
will  be  removed  from  the  fluid  by  the  external  loop  and  may  be  greater  than  or 
less  than  that  of  equation  (13)  depending  on  the  magnitude  of  various 
temperature  differences  and  changes.  Line  c'd'  represents  the  temperature 
profile  In  the  hot  column  after  heat  Is  removed  by  the  hot  end  external  loop. 

In  order  to  calculate  the  actual  cooling  output,  Qcold>  to  the  external 
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loop  and  the  system  coefficient  of  performance  (COP)  corrections  must  be  made 
to  heat  flow  calculated  for  the  Gadolinium  alone  in  equations  (10).  Based  on 
figure  8,  a  correction  for  the  cooling  output  was  made  where: 


Ql-2  *  Fcold 


where 


Fcold  =  1  ' 


atfc  ,  .  atgd. 


6-1  + 
— at: 


for  the  hot  end,  the  heat  load  to  be  removed  by  the  external  loop,  Q^ot,  was 
found  by  correcting  the  magnetic  heat  output  of  equation  (13)  by. 


^hot  ”  ^4 


where 


Fhot  =  1  + 


Thus  the  heat  load  to  the  external  hot  end  loop  can  be  greater  than  the 
magnetic  heat  load  of  equation  (13)  if  the  frictional  term  is  larger  than  the 
heat  lost  from  the  fluid  to  heating  the  Gadolinium. 


In  calculating  COP,  the  work  Input  to  the  loop  is  that  obtained  from 
equation  (13)  without  the  frictional  term.  This  value  does  not  include  any 
losses  In  generating  the  magnetic  field  which  would  Increase  the  actual  work 
input  to  the  system  necessary  to  generate  the  amount  of  cooling  given  by 
equation  (34).  Thus 
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t 


1 


'4-5  - 


JL  -o  _  . 

r)  ycold 

m 


where  is  the  efficiency  of  generating  and  transmitting  the  magnetic  field. 

CYCLE  RESULTS 

In  the  analysis,  it  was  desired  to  maintain  an  external  chilled  water 
loop  at  45°F  (7.2°C)  cold  end  with  an  external  hot  end  sink  temperature  of 
85°F  (29.4°C).  No  loss  effects  were  assumed  in  the  magnetic  cicuit  at  this 
time  or  in  the  external  cooling  loops.  Thus  only  the  performance  of  the 
Ericsson  cycle  was  analyzed.  The  effects  of  the  following  parameters  on  the 
cycle  performance  were  evaluated: 

o  Plate  Spacing,  Sp 

o  Channel  Length,  Lp 

o  Plate  Thickness,  tp 

o  Mixing  Temperature  Change,  ATmi-x 

o  External  heat  exchanger  temperature  differential, 

The  effects  of  channel  height,  HCH,  made  no  noticeable  change  in  cycle 
performance  for  the  thermal /hydraul ic  analysis  performed  and  was  not 
considered.  However,  channel  height  is  important  in  considering  the  magnetic 
circuit  since  it  represents  the  magnetic  flow  path  length  which  affects  magnet 
performance. 

Figures  9  to  13  show  the  effects  of  varing  the  above  parameter  for  a 
magnetic  field  strength  of  7-Tesla  while  figures  14  to  18  show  similar  results 
for  a  2-Tesla  field  strength.  Figure  9  through  18  are  plotted  for  cooling 
capacity  to  the  external  loop  and  coefficient  of  performance  versus  the 
Gadolinium  plate  channel  velocity,  Vp.  The  plate  channel  velocity  is  of 
importance  since  it  determines  whether  laminar  or  turbulent  flow  and  heat 
transfer  models  are  used.  The  figures  were  computer  plotted  using  a  spline 
fit  of  the  data  points.  Some  irregularities  in  the  curves  occurred  at  higher 
velocities  where  the  laminar  to  turbulent  change  occurs  between  the  plates  and 
the  prediction  models  change. 
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COOLING  CAPACITY  -  BTU/HR/LB  OF  GO  (SOLID  LINES) 


MAGNETO  THERMAL  HEAT  PUMP  PERFORMANCE  AT  7.0  TESLA 
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Figure  9a:  Effects  of  Gadolinium  Plate  Spacing  on  Performance  at  7-Tesla 
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Figure  9c:  Effects  of  Gadolinium  Plate  Spacing  on  Performance  at  7 -Tesla 
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Figure  10b:  Effects  of  Gadolinium  Channel  Length  on  Performance  at  7-Tesla 
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Figure  10c:  Effects  of  Gadolinium  Channel  Length  on  Performance  at  7 -Tesla 
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Figure  11c:  Effects  of  Gadolinium  Plate  Thickness  on  Performance  at  7-Tesla 
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Figure  12b:  Effects  of  Thermal  Mixing  on  Performance  at  7 -Tesla 


36 


DTNSRDC  PAS-82 -11 


COEFFICIENT  OF  PERFORMANCE  (DASHED  LINES) 


COOLING  CAPACITY  -  BTU/HR/LB  OF  GO  (SOLID  LINOS) 


MAGNETO  THERMAL  HEAT  PUMP  PERFORMANCE  AT  7.0  TESLA 

0*W*L  XISOT.r*.-  3.00  M.ATC  SPnCIM.lN.  -0.040  TCHf.KJT  no,r  •  tS.  N1X1NC  IDT,r-  2.0 
OMNI  LD44TH.IN.-  0.30  rt.ATC  PlJOUCM.  IN.-  0.010  TOT  COLD  CIC.r-  15.  N.C. OCLIB-T.f *  0.0 


3800.0 


3450. 0  + 


2100.0  + 


1750.0  + 


1400.0 


1050.0  + 


700.0  + 


350.0  A 


>10''  id1  10' 

GAOGLINIUM  PLATE  CHANNEL  VELOCITY  -  FT/SEC 


12.0 


r  U.O 


r  10. 0 


9.0 


r  8.0 


h  7.0 


6.0 


5.0 


P  4.0 


3.0 


r  2.0 


I-  1.0 


0.0 


Figure  12c:  Effects  of  Thermal  Mixing  on  Performance  at  7-Tesla 
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Figure  I3a:  Effects  of  External  Heat  Exchanger  Temperature  Difference 
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Figure  13c:  Effects  of  External  Heat  Exchanger  Temperature 
Difference  on  Performance  at  7-Tesla 
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Figure  14b:  Effects  of  Gadolinium  Plate  Spacing  on  Performance  at  2-Tesla 
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Figure  15b:  Effects  of  Gadolinium  Channel  Length  on  Performance  at  2-Tesla 
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Figure  16a:  Effects  of  Gadolinium  Plate  Thickness  on  Performance  at  2-Tesla 
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Figure  16b:  Effects  of  Gadolinium  Plate  Thickness  on  Performance  at  2-Tesla 
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Figure  17a:  Effects  of  Thermal  Mixing  on  Performance  at  2-Tesla 
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Figure  17b:  Effects  of  Thermal  Mixing  on  Performance  at  2-Tesla 
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Figure  17c:  Effects  of  Thermal  Mixing  on  Performance  at  2-Tesla 
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Figure  18a:  Effects  of  External  Heat  Exchanger  Temperature 
Difference  on  Performance  at  2-Tesla 
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Figure  18b:  Effects  of  External  Heat  Exchanger  Temperature 
Difference  on  Performance  at  2-Tesla 
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Figure  18c:  Effects  of  External  Heat  Exchanger  Temperature 
Difference  on  Performance  at  2-Tesla 
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COEFFICIENT  OF  PERFORMANCE  (DASHED  LINES) 


The  results  presented  in  figures  9  through  18  all  show  the  highest  COP's 
to  be  occurring  at  the  low  velocities  while  maximum  cooling  capacity  is 
occurring  at  higher  velocities.  In  general,  cooling  capacity  increased  with 
channel  velocity  until  frictional  effects  exceeded  cooling  gains  and  the 
cooling  capacity  dropped  off.  Increasing  the  cycle  rate  increases  the  cooling 
capacity  until  losses  associated  with  establishing  larger  finite 
temperature  differences  between  the  Gadolinium  and  fluid  to  effect  heat 
transfer  at  the  higher  cycle  rates  offsets  any  further  gains.  These  losses 
result  in  a  decrease  in  maximum  COP  with  increasing  cycle  rate.  The  maximum 
cooling  capacity  for  a  given  cycle  rate  essentially  coincides  with  the  maximum 
COP  for  that  cycle  rate. 

Comparing  the  7  and  2-Tesla  systems  for  similar  geometries  and 
conditions,  the  7-Tesla  system  provide  approximately  4  to  6  times  the  cooling 
capacity  of  the  2-Tesla  system.  The  7-Tesla  system  yields  higher  COP's  than 
the  similar  2-Tesla  system  operating  at  the  same  plate  channel  velocity.  The 
7-Tesla  system  can  operate  at  higher  plate  channel  velocities  and  shows 
increasing  cooling  capacity  up  to  approximately  twice  the  cycle  rate  of  the 
similar  2-Tesla  system. 

Figures  9  and  14  show  the  effects  of  plate  spacing,  Sp,  on  the  system 
performace  at  7  and  2-Tesla,  respectively.  Increasing  plate  spacing  increases 
the  COP  and  permits  operation  over  a  wider  range  of  plate  channel  velocities 
with  maximum  cooling  capacity  being  shifted  to  higher  velocities.  Maximum 
cooling  capacity  increased  slightly  for  the  7-Tesla  case  but  started  to 
decline  some  for  the  2-Tesla  case.  Thus  increasing  plate  spacing  improves 
performance  of  the  cycle.  However  spacing  the  plates  further  apart  can  have 
an  adverse  affect  on  the  magnet  performance  because  of  the  increased  void 
facing  the  magnet  and  its  effect  on  magnet  performance  which  is  not  considered 
in  the  figures. 

Figures  10  and  15  show  the  effects  of  the  Gadolinium  channel  length,  Lp, 
on  the  performance  for  a  7  and  2-Tesla  system,  respectively.  Increasing 
channel  length  causes  a  significant  reduction  in  both  cooling  capacity  and  in 
COP  due  to  the  added  frictional  losses  in  the  longer  channel.  Cycle  rate  at 
the  maximum  cooling  capacity  also  decreased  noticeably  with  increased  channel 
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length.  Thus  the  shorter  the  Gadolinium  channel  length  the  better  the  cycle 
performance. 

Figures  11  and  16  show  the  effect  of  the  Gadolinitmi  plate  thickness  on 
performance  for  a  7  and  2-Tesla  system,  respectively.  Reducing  plate 
thickness,  significantly  increased  cooling  capacity  and  enabled  maximum 
cooling  capacity  to  occur  at  higher  cycle  rates.  Maximum  COP's  remained  about 
the  same  for  any  given  plate  channel  velocity  but  occurred  for  slightly  higher 
cycle  rates  at  reduced  plate  thicknesses.  The  thinner  plates  improve  heat 
transfer  because  of  less  thermal  resistance  and  provide  increased  surface  area 
per  pound  of  Gadolinium  which  lowers  finite  temperature  differences  between 
the  fluid  and  ,:late  necessary  for  heat  transfer,  but  increases  frictional 
losses  due  to  the  increased  area.  However,  thinner  plate  thicknesses  are 
desirable  for  improving  cycle  performance  neglecting  the  magnet  performance 
which  could  be  adversely  affected  by  facing  a  larger  void  with  the  thinner 
plates. 

Figure  12  and  17  show  the  effect  of  thermal  mixing  on  the  perforance  for 
a  7  and  2-Tesla  system  respectively.  Thermal  mixing  results  from  such  things 
as  boundary  layer  formation  in  the  column  causing  warmer  fluid  to  be  moved 
into  a  region  of  cooler  fluid  and  vise  versa.  Thermal  mixing  was  expressed  as 
a  temperature  difference  along  the  cold  and  the  hot  column  caused  by  the 
thermal  mixing  and  was  assumed  to  be  the  same  in  both  end  columns.  Increased 
mixing  temperature  decreases  cooling  capacity  slightly  and  reduces  COP 
noticeably.  Increased  mixing  temperature  results  in  more  of  the  cooling  load 
being  required  to  be  used  in  regeneration  to  cool  the  Gadolinium  plates  as  the 
cold  column  passes  through  the  plates  (process  6  to  1)  prior  to  the  reverse 
trip  (process  1  to  2)  where  the  magnetic  field  is  dropped  to  obtain  the 
cooling  load.  Thus,  less  of  the  cooling  load  is  available  for  the  external 
cooling  loop. 

Figure  13  and  18  show  the  effect  of  the  heat  transfer  temperature 
difference  between  the  external  loop  and  the  cycle  fluid  for  a  7  and  2-Tesla 
system,  respectively.  Increasing  heat  transfer  temperature  difference  reduces 
both  cooling  capacity  and  COP.  This  results  from  the  cycle  operating  over  a 
larger  temperature  range,  thus  more  heat  must  be  transferred  in  the  same 
amount  of  time.  Thus,  larger  heat  transfer  temperature  differences  are 
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required  between  the  fluid  and  the  Gadolinium  which  result  in  the  performance 
decrea se. 


SYSTEM  SIZING 

An  important  consideration  in  selecting  a  system  is  its  physical  size. 

To  avoid  gravitational  mixing  due  to  differences  in  fluid  density,  the  column 
would  operate  vertically  with  the  hot  column  on  top  and  the  more  dense  fluid 
of  the  cold  column  on  the  bottom.  The  system  size  is  primarily  a  function  of 
plate  geometry  and  plate  channel  velocity.  It  is  assumed  that  the  frontal 
area  of  the  plates  is  the  same  as  the  cross-sectional  area  of  the  column. 

Thus  the  flow  through  the  plates  is  at  a  higher  velocity  than  the  column 
velocity  due  to  the  reduced  free  flow  area  caused  by  plate  blockage.  The 
ratio  of  free  flow  area  through  the  plates  to  the  column  cross-sectional  area 
is  expressed  by 


R 


A 


t  +  S. 


(39) 


and  the  cross  sectional  area  of  the  column,  AqL  is  found  by 


GO 


12 

n-^rr 


P 


(40) 


figure  19  shows  the  effects  of  the  free  flow  area  ratio,  RAj  and  plate  channel 
length  on  the  column  cross-sectional  area.  Very  short  plate  length  and  high 
free  flow  area,  which  show  good  performance  trends,  are  high  in  cross- 
sectional  area  and  thus  volume. 


The  fluid  column  length  represents  about  half  the  system  length  since  the 
column  must  be  moved  completely  through  the  plates  and  back  again  to  complete 
one  cycle.  System  length  was  calculated  from 
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GADOLINIUM  PLATE  LENGTH  -  INCHES 


Figure  19:  Column  Cross  Sectional  Area 


60 


DTNSRDC  PAS-82 -11 


(41) 


L 


s 


60  RA  VP 

nrpm 


Figure  20  shows  a  plot  of  system  length  as  a  function  of  plate  channel 
velocity,  cycle  rate,  and  free  flow  area  ratio  assuming  plate  length  equals 
zero.  Increasing  cycle  rate  decreasing  free  flow  area,  and  reducing  plate 
channel  velocity  will  reduce  system  length.  Multiplying  equations  (40)  and 
(41)  would  give  an  idea  of  system  volume  per  pound  of  Gadolinium  not  including 
the  magnetic  portion  o*  the  system.  System  volumes  will  be  larger  than 
conventional  systems  to  Improve  performance. 

DISCUSSION  OF  RESULTS 

Present  room  temperature  refrigeration  systems  operate  with  best  COP's  on 
the  order  of  3  to  4.  For  Gadolinium  magnetic  heat  pumps  to  be  of  interest 
COP's  of  at  least  5  or  greater  would  be  necessary  to  provide  any  improvement 
over  present  state-of-the-art  systems.  Reviewing  figures  9  to  18  show  that 
COP's  of  5  or  greater  occur  at  plate  channel  velocities  of  2.5  FPS  or  less  for 
a  7-Tesla  system  and  at  1.5  FPS  or  less  for  a  2-Tesla  system.  At  these 
conditions  the  7 -Tesla  system  showed  a  maximum  cycle  rate  of  12  cycles  per 
minute  and  a  maximum  cooling  capacity  of  about  690  BTU  per  hour  per  pound  of 
Gadolinium.  For  the  2-Tesla  system,  maximum  cycle  rate  was  4  cycles  per 
minute  with  a  maximum  cooling  capacity  of  about  95  BTU  per  hour  per  pound  of 
Gadolinium.  At  higher  plate  channel  velocities,  maximum  cooling  capacities  of 
as  much  as  475  and  2300  BTU  per  hour  per  pound  of  Gadolinium  were  shown  for  a 
2  and  7-Tesla  system,  respectively ,  but  COP  was  less  than  two. 

Figures  9  to  18  presented  the  effects  of  the  various  parameters  on  the 
system.  Based  on  this  information  a  more  optimum  system  was  evalauted  using 
what  were  felt  to  be  reasonable  geometries  and  losses.  The  system  shown  in 
Table  2  was  evaluated  and  presented  in  figures  21  and  22  for  a  7  and  2-Tesla 
system  respectively. 
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Table  2  -  Optimum  System  Parameters 
Field  Strength  7  Tesla  2  Tesla 

Plate  Spacing,  Sp  .060  Inches  .060  Inches 

Plate  Thickness, tp  .005  Inches  .005  Inches 

Plate  Length,  Lp  .250  Inches  .250  Inches 

Mixing  Temperature,  ATmi-x  0.1  °F  0.1  °F 

Heat  Exchanger,  atR£  1.0  °F  0.5  °F 

Using  the  more  optimum  plate  geometry  even  with  some  mixing  temperature 
effects,  system  performance  was  improved.  External  heat  exchanger  temperature 
difference  was  assumed  higher  for  a  7-Tesla  system  since  it  can  operate  at 
higher  cycle  rate  which  requires  the  larger  value.  Figures  21  and  22  show  the 
performance  of  the  7  and  2-Tesla  systems  respectively  of  Table  2.  Table  3 
tabulates  the  optimum  performance  as  well  as  its  column  size  for  the  system 
based  on  a  COP  of  5.  Maximum  cooling  capacities  at  the  COP  of  5  turned  out  to 
be  1200  BTU/HR/LB  of  GD  at  7-Tesla  and  130  BTU/HR/LB  of  GD  at  2-Tesla. 


Table  3  -  Optimum  System  Performance  &  Size 


Field  Strength 

7  Tesla 

2  Tesla 

COP 

5 

5 

Plate  Channel  Velocity,  Vp 

3  FPS 

2  FPS 

Cycle  Rate,  NRpM 

22  cpm 

6  cpm 

Cooling  Capacity,  QCo1d 

1200  BTU/HR/LB  of  GD 

130  BTU/HR/LB  of  GD 

System  Length,  Ls 

7.6  ft 

18.5  ft 

Column  Cross-Sectional 

Area  (ACL)  per  ton  of 
refrigeration 

12.7  ft  2/ton 

117.4  ft  2/ton 

Column  Volume  per  ton  of 
refrigeration 

96.5  ft  ^/ton 

2172  ft  3/ton 

Higher  cooling  capacities  are  possible  but  only  at  COP's  less  than  or 
comparable  to  present  state-of-the-art  systems  and  thus  no  performance 
Incentive  to  consider  magnetic  heat  pump  application. 
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COOLING  CAPACITY  -  BTU/HR/LB  OF  GO  (SOLID  LINES  I 


MAGNETO  THERMAL  HEAT  PUMP  PERFORMANCE  AT  7.0  TESLA 
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Figure  21:  Improved  Performance  of  a  7-Tesla  System 
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COEFFICIENT  OF  PERFORMANCE  (DASHED  LINES) 


COOLING  CAPACITY  -  BTU/HR/LB  OF  GO  (SOLID  LINES) 
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Figure  22:  Improved  Performance  of  a  2-Tesla  System 
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COEFFICIENT  OF  PERFORMANCE  (DASHED  LINES) 


The  1200  and  130  BTU/HR/LB  of  GO  for  the  above  system  (COP  =  5)  is  the 
maximum  to  be  expected  with  any  significant  improvement  in  COP.  Mixing 
temperature  loss  was  assumed  to  be  0.1°F  in  the  above  system.  Low  mixing 
temperature  loss  may  only  be  possible  if  velocity  profile  formation  in  the 
column  is  retarded.  This  could  require  some  arrangement  such  as  walls  that 
move  with  the  column  to  avoid  velocity  profile  formation  across  the  column. 

The  above  system  will  require  10-pounds  of  Gadolinium  to  obtain  one- ton  of 
refrigeration  (12,000  BTU/HR)  at  7-Tesla  and  92.3  pounds  of  Gadolinium  at  2- 
Tesl  a. 

Table  3  also  includes  an  idea  of  the  column  size  obtained  from  use  of 
equations  39  to  41.  Column  volumes  of  96.5  and  2,171  ft3/ton  of  refrigeration 
were  obtained  for  a  7  and  2-Tesla  system  respectively.  Column  lengths  were 
7.6  and  18.5  ft  for  the  7  and  2-Tesla  system  respectively.  These  systems  were 
sized  for  optimum  performance  rather  than  minimum  size.  Size  could  be  reduced 
some  at  the  expense  of  performance  but  would  be  considerably  higher  than  the 
estimated  4-8  ft3/ton  for  a  freon  type  air  conditioning  system.  The  above 
column  volumes  do  not  include  the  volume  of  the  magnetic  system.  The  magnetic 
heat  pump  system  will  be  considerably  more  voluminous  than  the  more 
conventional  state-of-the-art  systems.  Unless  the  field  strengths  could  be 
pushed  considerably  beyond  the  7  Tesla  limit  considered,  magnetic  heat  pump 
column  size  can  not  be  significantly  reduced  to  compete  with  conventional 
systems. 

In  this  analysis,  the  magnetic  circuit  efficiency,  r?  ,  was  taken  as 
unity.  Any  losses  in  the  magnetic  circuit  are  not  reflected  in  the  system 
OOP's  that  are  presented.  Generating  the  magnetic  field  represents  the  work 
input  to  the  system  which  is  part  of  the  basis  on  which  COP  is  determined.  In 
equation  (38)  for  C0P,fQ4_5-Qp^  ^represents  the  energy  drawn  from  the  magnet 
by  the  Gadolinium  plates.  Any  losses  in  generating  and  transmitting  this 
energy  to  the  plates  is  included  1^  the  magnetic  circuit  efficiency  and  makes 
the  work  input  to  the  magnet  equal  I Q^_g  -Qp4  5) /’’m  •  ^9ure  23  shows  the 
effect  of  magnet  circuit  efficiency  on  COP  calculated  from  equation  (38) 
compared  to  an  ^  =1  ideal  case.  For  the  COP  of  5  case,  a  magnetic  circuit 
efficiency  of  0.9  would  result  in  an  actual  COP  of  only  3.0  which  barely 
approaches  state-of-the-art  systems.  Magnetic  circuit  efficiencies  of  greater 
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COEFFICIENT  OF  PERFORMANCE  (ACTUAL) 


Figure  23:  Effects  of  Magnetic  Circuit  Efficiency  on  COP 
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than  .95  appear  to  be  necessary  to  assure  high  COP  performance. 

The  magnetic  heat  pump  has  been  analyzed  starting  with  a  series  of  ideal 
processes  and  modifying  these  for  the  major  losses.  Friction  losses  could  be 
calculated  with  a  reasonable  degree  of  certainty.  However,  mixing  losses  in 
the  column  were  difficult  to  evaluate  and  were  estimated  over  a  reasonable 
range  to  show  its  effects.  Actual  calculated  estimates  for  mixing  losses 
would  be  required  for  any  further  serious  consideration  of  the  cycle.  The 
cycle  was  evaluated  with  magnetic  field  strengths  up  to  7 -Tesla  which  is  the 
present  limit  of  superconducting  magnets.  Improvements  in  field  strength 
above  this  limit  should  impact  favorably  on  this  magnetic  heat  pump  cycle  by 
Improved  performance  per  pound  of  Gadolinium  and  reduced  columrn  size.  The 
external  heat  exchanger  loop  which  removes  the  cycle  cooling  load  would 
ideally  operate  over  a  very  small  heat  transfer  temperature  difference  to 
improve  cycle  performance  and  the  impact  of  the  actual  heat  exchanger  design 
on  the  system  must  be  considered  in  any  further  consideration  of  cycle. 


CONCLUSIONS 

Based  upon  the  idealized  thermal  analysis  of  the  magnetic  heat  pump  the 
overall  conclusion  is  that  the  particular  system  studied  is  not  attractive  for 
shipboard  use  due  to  its  large  size  and  unspectacular  perfomance.  An  improved 
analysis  and  consideration  of  the  magnetic  and  electrical  portion  of  the 
system  will  most  likely  show  that  the  performance,  weight  and  volume  are 
inferior  to  conventional  freon  heat  pumps. 

The  following  specific  conclusions  should  be  noted: 

o  A  7-Tesla  system  would  provide  approximately  4  to  6  times  the 
cooling  capacity  of  a  2-Tesla  system  for  similar  conditions  and  with 
COP's  as  much  as  30  to  50  percent  higher. 

o  The  ideal  plate  geometry  would  minimize  plate  thickness  and  plate 
length  and  increase  plate  spacing  to  provide  best  cycle  performance. 

o  At  a  cycle  COP  of  5,  a  7 -Tesla  system  could  provide  a  cooling 
capacity  of  1200  BTU  per  hour  per  pound  of  Gadolinium  while  a  2 
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Tesla  system  would  operate  around  130  BTU  per  hour  per  pound  of 
Gadolinium.  Higher  cooling  capacities  are  possible  but  at 
unfavorable  OOP's  while  higher  OOP's  are  possible  but  at  reduced 
cooling  capacities. 

o  Magnet  circuit  efficiency  must  be  high  to  obtain  high  OOP's. 

Magnet  circuit  efficiency  of  80  percent  or  less  would  result  in 
actual  operating  OOP's  that  are  lower  than  best  conventional  systems 
that  range  around  COP  equal  3  to  4.  Magnet  circuit  efficiencies  of 
95  percent  or  better  are  needed  to  consider  magnet  heat  pumps  for 
improved  performance. 

o  The  magnetic  heat  pump  investigated  would  be  considerably  larger 
than  a  conventional  state-of-the-art  cooling  system  by  as  much  as  an 
order  of  magnitude. 

o  Increasing  magnetic  field  strength  can  decrease  the  system  size 
significantly  by  increasing  the  amount  of  cooling  capacity,  possible 
per  pound  of  Gadolinium.  However  magnetic  field  strengths  much 
greater  than  the  present  superconducting  state-of-the-art  limit  of 
7-Tesla  must  be  possible  before  the  magnetic  heat  pump  investigated 
could  compete  in  size  with  conventional  systems. 


69 


DTNSRDC  PAS-82 -11 


REFERENCES 


1.  Giauque,  W.  F.  and  0.  P.  MacDougall,  Phys.  Rev.  43,  768,  1933 

2.  Steyert,  W.A.,  "Sterling-Cycle  Rotating  Magnetic  Refrigerators  and  Heat 
Engines  for  use  near  Room  Temperature",  J.  Appl .  Phys.  Vol .  49  No.  3, 

March  1978 

3.  Brown,  G.V.,  "Magnetic  Heat  Pumping  Near  Room  Temperature",  J.  Appl.  Phys. 
Vol.  47  No.  8,  August  1976 

4.  Norris,  R.H.  and  0.0.  Streid,  "Laminar-Flow  Heat-Transfer  Coefficients  for 
Ducts",  Transactions  of  the  A.S.M.E,  August  1940 

5.  Kreith,  F.,  Principles  of  Heat  Transfer,  International  Textbook  Co.,  April 
1963 

6.  Brown  G.V.  and  S.S.  Papell ,  "Regeneration  Tests  of  a  Room  Temperature 
Magnetic  Refrigerator  and  Heat  Pump,"  Preprint  Submitted  to  Applied 
Physics  Letters,  March  1978 


70 


DTNSRDC  PAS-82 -11 


Appendix  A 


Magnetic  Heat  Pump  Performance  Computer  Program 
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o  o  o 


PO0r,D4M  HAGHEAT  (INPUT,  OUTPUT,  T APES = T NPU T , T AO E6=0UTPU T ) 

THIS  PROGRAM  =.VALUATcS  THE  PERFORMANCE  OP  AN  ERICSSON  C  YCL  c 
MAGNETIC  HEAT  PUMP  FOP  ROO“  TEMPERATURE  APPLICATIONS  USING 
GADOLINIUM, 

COMMON  OCX  (50 ,50  )  ,  VELX  (SO  ,6  0)  ,.NRX,NV(50)  ,  COPX  (  5  0 , 50  ) 

COMMON  WC,LC,TC,TP,TFH,Tei-;,DTM,QTW,OCU 
COMMON  N , NN 

RrAL  LC,KW,KF,MUF,NUF,LrN,PR 
INTEGER  PHI 

DIMENSION  PHI (32) , NUF (3?) ,  TCX (10) ,  RHOX (101 , CF  X  (  10  ) , XKF (lOt , XMU (10) 
•  *  T ( 3 ) , FQ ( 1 1 ) , V(30) 

OAT A ( (NUF ( I) ,1=1, 32) = .5, 1. , 1.5 ,2., 2. 48, 2.35,3.42, 3.85,4 . 23, 4.58, 
*5. 13, 5. 89, 6.  67, 7. 16, 7. 53,o. "A, 8. 67, 9.  5*, 19. 4,  11.6,13.1,14.4, 16.2, 
*19.  5,23.9,26.  9,29.3,3  3.  5, 39. 9,50.2,57. 2,62.8) 

DATA ( <  PHI (J)*J=1»32)  =  1, 2,3, 4, 5, 6, 7, 3,  9,  10,12,  15,20, 25,33, AO, 60, °0, 
*130, 203, 300, 400, 600, 100), 2000,3000, 40 00, 60 00, 13000,20000, 30 000, 
*40000) 

DATA  < (TFX(I) ,1*1, 9) *3 2., 40. , 50 . , 60 . , 7 0 . , 8 0 . ,9 0. ,  10  0.,  15  0.) 

DATA ( (RmOX (I)  ,1  =  1, 9) =62. 42, 82. 42, 62. 42, 6?. 3 8, 62. 31, 62. 27, 62. 11  , PI. 
*88,61.20) 

QATA( (CFX(I) ,I  =  1,9)  =  1.01, 1.  )3.1.0Q» .9  99,  .998,  .998, . 397, .998, 1  .0  0) 
OATA ((XKF (I), 1=1,9)=. 3 19, .325, .332, .3 40, .347, .353, .359, .3 64, .784) 
DATA ((XMU(I) ,1=1,91=1.2, 1.C4,. 88, .76, .658, .573, .514, .458, .292) 
OATA((FQ(J),J=l,ll)sl.,2.t4.,3.,12.,18.,24.,30.,40.,50.,60.) 

DATA  (  ( v/  (  K )  ,K=1,2Q  )  =  .  5, 1.  ,2,  *  3.  ,4., 5.  ,6.  ,7.  ,3.  ,9.5,11.  ,12,5,14.  ,  15. 
15, 17., 18. 5, 20. ,21. 5 ,22.5,23.5) 

N=  0 


NN=  1 

1  REAP (5 ,100) WC,TC,LC,TO,9TM,9Cn,9TW 
100  FORMAT (3F13. 4) 

N*>X=  7 


IF(WC.EQ.O.Q)  GO  TO  999 
DO  60  J=l,ll 
FR"0=FQ ( J) 

WRITE (6,101)  WC»LC* TC»  t° 

101  FORMAT  (1H1/2X,  “CHANNEL  WIDTH  =",F6.3,”  IN" , 5X ," CH AN NEL  LENGTH  =",F 
*7.3,"  IN", 5X, "PLATE  SPACING  =",F6.3,"  I N~  ,5 X , "PL A TE  THICKNESS  =",f 
*6.3,"  IN"/) 

TFH=  , 


TFC=45. 

TFI= (TFH^TFC) *,5 
T ( 1 )=TFC-OTW 
T (2) =TFh»OTW 

T(?I=TFX 

TIRP=1  . 

9HU=9CU 
PCL  =  0. 

CALL  f  DPROH ( 1 ,8CU, T (1) ,OTCU,STCUI,BTCU) 
CALL  GOPROPCl,BCL,T(l) ,9TCL ,STCLI,BTCL) 
CALL  GOPROP (1 ,3MU,T (2 ) ,9  THU, STHUI, BTHU) 
STHLI=STHUI*STCLI-STCUI 
S TO T= (STHUI -STCUI I  /  <  T  <  2 )  -T  ( 1  >  > 

AFx«r*TC/ 144. 

Af»mc*LC/144,*2. 

0t=4.*LC/l2.*AF/AS 
0H*2.*WC*TC/ (TC*WCI /I 2 . 

RHOW=  490 • 7 


^OV 
r  &°C  >-,o(V 
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I 


i  j l  *"  :±?sir 


i 


l 

! 

I 


< 


I 


i- 


KW=  5  .  * 


C°GO=0.07l 


NPASS=1728./RH0W/ <WC*LC*TP)  ♦  .1 

G0MCLE=157 .26 

70=459.7 

HR  I TE<6, IQ?)  TFH,TFC,P9?0fNoasS,0TM, OTM.gru 
132  FPPMA  T  <2X,**  TEMPERATURE  MOT  ENO  =  **,F6.lt-  f  TEMPERATURE  COLO 

*0  =  **,F6.i,*‘  F  CYClES/MINUTE  =“,F5  .t  ,5X,**N0.0F  CHANNELS/L9  =**,I 

*5//2X,**H£AT  EXCHANGER  TEMPERATURE  DIFFERENCE.  =“,F5.1,“  F-,5x,“MIxi 

*ng  temperature  change  =**,F6.i,-  f**,tx,**  max.fielo  =-,P4.i,**  tesla** 
*//) 

W0IT^(6, 198) 

WPITF16* 19P) 

198  FORMAT  (IX,**  VEL  REV. NO.  HAF  TIME  OEL-T  PHI  F 

*HL  U  UA  QP  QF  OT-CHR  t  I  HE  LENGTH*** 

199  F  OP  MAT  ( IX  ,  ’*  (FPS)  (SEC!  (IN 

*.H20)  (9/HF2R) (a/H-R )  (9/H)  (9/H)  (F)  (SEC)  (FT)-) 

NVX=0 

CO  60  K=l,20 
CCX=n. 

OHX=0. 
dcf  =  a . 

PWF=A. 

OFF =  0  . 

qfh=g. 

CFFC - Q  » 

HFFHiQ  . 

OGDC  =  0  . 

OGDH=Q  . 

NCT  =  0 


OTPF  G*0  . 

T I RP  =  1  • 

3  VCL=V(K) 

VELCH=VEL*TC/  (TOTP) 

CO  4P  L*l.  3 
no  5  1=2,9 

IF(t(L).LT.TFX(I) )  GO  TO  $r 

6  CONTINUE 

30  FP0= (T(L) -TFX(I-l) )/ (TFX ( I) -TFX(I-1) ) 
°HOF  =  RHOX <1-1 ) * (1 ,-FRC)  *9H0X  (I)*FBC 
CF=CFX(I-1) * (1.-F9C) ♦CFX (I) *FPC 
KF=XKF(I-1) *  11 ,-FBC) ♦XKF(I)*F9C 
MUFaXMUl  I-l)  *  (l.-FBC)  ♦  XM(J  (  I )  *F9C 
MUF=  MUF*  Q. 001 
PP=ff*MuF/<F*3630 . 

GF=RHOF*VEL*3600. 

ret=phof*vel*oh/muf 

IF ( PET .GT.3000.)  GO  TO  7 
PH=cF*GF*OE*0f/ (<F*LC/12. ) 
F=96./(R£Y*(1.*TC/WC>  **2) 

GO  TO  8 

7  Fs0.184/RE v**o,2 

8  HL=LC/0H*VEL*VEL/2.*F/32.174 
IF (REY.LT .3000.)  GO  TO  21 

HAF  =0. 023*KF/OH*RE Y* *  3 . 8 *PR**0. *3 
GO  TO  22 
21  CONTINUE 


,o$< 
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no  10  1  =  2,  *2 
IFtPH-PHUIl)  IS, 20,  10 
10  CONTINUE 

15  FAC=  (3H-PHI  (I-m  /  (PHI  m-PMl  (1-1)  1 
FNU=NUF(I-1I  ♦FAC*  (NUF  (I)  -Nl)P(I-ll) 

GO  to  25 
20  FNU=NUF(II 
25  HAF=FNU*KF/OE 

22  TIM^sUC/ ( 12.*VEL) 

TlMH  =  TIHS/3600  . 

U=l./ (1./HAF+TP/24./KW) 

UAS=U* AS 

IFtTIPR.EQ.l.)  GO  TO  23 
IF (L-2)  26,27,25 

23  0TWC= (OTM+OCF) /2. 

TCU=T(i) -OTWC 
TCUW=T (1 ) ♦OTWC 
TCLMsTCUM-OFC 

TCL  =  TCLW-QCX 

IF(OGOC.GT.OFC)  GO  TO  55 

CALL  GOPROP ( 1 «  9CU, TC U , O TCU, S TCU , 9  TCI) ) 

CALL  GOPROp (1 , 3CL , T  CL , OTCL , STCL , 0TCL ) 

IF (STCL,LT  .STCU)  GO  TO  50 
OTC=nTCU-OTCL 
STC=STCU-STCL 
OTWH= (OTN-OHF) /2. 

Thl=T (2) +  PTWH 
THLH=T<2)-0TWH 

thuw=  thlw+ofh 
thu=thuw+ohx 

CALL  GOPROP (1 ,dHU,THU,0TMU,STHU,9TrtU> 

STHL  =  STHLI*STOT* (THL-T(2)  ) 

IP (STHL.LT.STHU)  GO  TO  50 

CALL  GOPROP (2.STHL, THL, OTHL , STHL  *9  THL > 

OTH=OTHU-PTHL 
STH=STHL -STHU 

OPG=( (THU-TCUH-uTREGI ♦ (THLW- JTREG-TCL > ) *CPGC*PREQ*60. 
UR=0»G 

QG  = ( STCL-STCU) /GOHOLE*  C  <TCL  +  TCU)/2.  +459.7) *FPEQ*60 . 
QH= (STHL-STHU) /GONOLE* (tTHL  +THU) 7  2.  *459. 7 ) *F RE Q*6 0 . 
QTOT  =  q(?«.<  jO+QH!*2. 

TICC=QC/ QTOT 
T IRR=QR/ QTOT 
tihh=qh/qtot 
GO  TO  T 
’6  QO=OC 

TTP=TICC/FREQ*60. 

0TR«-G=QC/UAS/TICC/NPASS 

OCX=OTREG 

LEN=VFLCH*TIM 

clcn=len 

CW= ISTCL-STCU) /GOMOLE* (t (1) +459.7) /OTC 
OFLT=FXP  <-HAF*TIMH/ <CW*R HOW* TP/24. ) ) 

QF=HL*RHOF»AS*VEL  *4. 6  295*TICC/12. *NPASS 
QGF=QF 

HCNT=GF*AF*TICC*CF*NPASS 

OFFC=OF/HCNT 


'  ■  C  does  bo) 
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QFC  =  (QC-QF) /HCNT 
IF<NCT  .LT.il  GO  TO  40 
WPI TF<6*  103) 

103  FOP MAT (1H  ) 

AP£P= (OFC-DFFC-OGOC) /OFO 
IF (APER.L£ • 0. )  GO  TO  50 
QP= (QC-QF) *APpP 
GO  TO  20 

27  QP=OH 

TIM=TIHH/FP£Q*60. 

OTRFG=JH/UAG/TIHH/NPASS 
OhX  =  OTREG 
LEN= VELCH*TIH 
HLEN=LEN 

CW= (STHL-STHU)  /GO  HOLE  *  (T  (  2 )  *-4  50.  71/0  TH 

0ELT=£XP<-4AF*TIHH/(CW*pH0W*Td/ji,.)) 

QF=hL*RH0F*AS*V£L*4.6293*TI 4H/12.*NPASS 
GHF  =  OF 

HCN  T=GF*AF*TI HH*CF*NP A  SO 

offh=qf/hcnt 

0FH= (QH7QF) /HCNT 
IFtNCT.LT.T)  GO  TO  40 
AP£rl=  (OFhOFFh-OGOHI  /CFH 
GO  TO  29 

28  QRaQPG 
TIH=TTRR/FREO*60. 

LfN=VcLCH*TIM/2. 

RLFHsLEN 

CTREG=QR/UAS/TIPR/NPASS 

CW= ( STHU-STCO) /GOMOLE*  CTPT*45o. 7) / ( TFH-TFC) 
°EL  TaFXP (-HAFMIMH/ (CH*°H0W*To/24.  J  I 
qpc-hl*RH0F*AS*V£L*4. 62P3*T Tor/12. *NPASS 
QFsORF 

HCNT=GF*AF*TIR9*CF*NPASS 

OCF=ORF/HCNT 

OHF=OCF 


OLOC=CPGO*FREQ*60.* (2 • *QTWC ♦CTREG) /{ &F* AF» T ICC *CF *N PASS • 
OGOH=CPGO*FPc.Q*60.  *  (  2  .  *  OT  WP  *-DTRE  G)  /  <  GF  *AF*T  I HH»CF  *NPASS  ) 

IF ( NC  T  ,L  T , 3 )  GO  TO  40 

29  MpIT£ (6.200)  VEL.REY, HAF,TIHSfOELT,PH,FfHL,U.UAS.  QR.QF, 

•OTREG, TIM, LrN 

2  L  0  FOPMAT(F6.2,2F9.2,F6.3,F7.t,F8.2,F8 .a,F7.3,F9.2,F8.3,F7.1,F7.2,F7. 

•3  ,F7 . 2  .Ftt  ,  3 )  Cj 

4  0  CONTINUE  XV 


48 


IF(NCT.LT.3)  GO  TO  48 
PETAaJTFC 7459. 71/ (TFH-TFC) 

COP= (OC-QCF) *ApFR/ ( QH- ( QC-CCF ) * A  PER ) 
PATTO=COP/8ETA 

TLEN  =  CLEN7HLEN7RLEN*LC/12. * (l.-TC/ (TC  7TP) ) 
TILEN=TLEN*12 . 

TLPATsTILEN/LC 
go  TO  49 
NCT=NCT7l 


cV 

<? 


\0  V 
.<$? 


\e 


T IRPs 1 . 

GO  TO  3 

49  WPITE(6,300)  BTHL,OTC,OTH,PCU,dETA,COP, RATIO 
300  FORMA  T  <  /5*,”8HL  =“ , F 6 . 3 , 5X , "OTC  =~ , F7 .2 f 5X ,”OTH  =" ,F7 . 2 , 5X  ,”9CU  = 
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*“,F6. 3 ,5X,"BETA  =" , F7 . 2 , 5X , "COP  =“ , F7 . 3 , 5 X, "RAT  10 < C OP/BET  A)  =",F7.4 
*) 

W°ITE(6,301)  STCU»STCL»STC,CTHU»STHL»STH,TLRAT 
301  FORMATJ  5X  *"STCU  =  " ,  F7 . 3 , 4  X,  "STCL  =", F7. 3 , 4X, "STC  =  ",  F7.  3, 4X,  "STHIJ 
*  =“,F7.3,4X,"STHL  =“,F7.3,4X,“STH  =  " ,  F7 . 3 , 5 X,  "L  ENGTH  RA TI 0=" , F8 . 2 ) 
WRITE  (6t  302)  TCL,  TCUt  THL  ,  THU,  TLEN,  TIL^ 

30  2  FORMAT (  5X,“TC (LB >  =  “ , F7 . 2 , 5 X ,"TC (UB) =“, F 7, 2 , 5 X , "T H ( LB) =  “ , F7 . 2 , 5X  , 
♦"TH<U8)=“,F7.2,5X,“T0TAL  LENGTH  =",F7.3,"  FT  (",F7.2,"  IN)") 

WRITE  (6,  104)  TCU,TCLW,TCL,CCF,OFC, OFFC , 0G0C,0 TWC, APER 
WRITE (6, 104)  THU, THUW, THL , OHF, OFH, DFFH , DGOH , 0 TWH, APEH 
104  FORMAT (9F14. 4) 

IFrAPER.LT. 0.)  GO  TO  60 
IFJQCF.GT.QC)  GO  TO  60 
IF ( NVX.EQ. 0)  NRX=NRX*1 
NVX=NVX>1 
VELX  <NRX,NVX)=VEL 
OCX  (NRX, NVX) s (QC-QCF) *ApER 
COPX(NRX,NVX)=COP 
NV(NRX)=NVX 
FO  CONTINUE 
60  CONTINUE 
N=N  +1 

CALL  GRNBRO 
GO  TO  1 
999  STOP 
RNO 


TlC  does 
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SUP  °CUT I Nf  G0P90P INCASE, x ,Y, rr,ST,9T) 

THIS  SUBROUTINE  CON  T  4 1  OS  Physical  ppoocpTI^S  OF  GADOLINIUM. 

01  HENS  ION  Tk ( 14) ,B(5) ,DT<(5, 14) ,S(5, 14) , DTL (141 ,SIL (14) 

CASe  1=  X  ( T  £  SL  A )  Y(TE*°)  CASF  2=  X(CNTROPY)  v(TFMP) 

CAT  A  ( (TK  (I)  ,1  =  1,  141  =  263.  ,  261.  ,  273.  ,  27  8.  ,28*.,  2  88.  ,293.  ,298.  ,  103.  , 
*709. ,313. ,318. ,323., 328. I 
DATA ((9(1), 1  =  1,5) =7., 5.,  3.  ,t.,Q.) 

CAT  A (( (DTK (J, I), 1  =  1,14) , J=l,5)» 


y 

7.40, 

8.00, 

O.70, 

9.65, 

10.90, 

12.31, 

1 3.95, 

13. 72, 

13.28, 

12.51, 

11.63 

« 

ir.95. 

10.20, 

9.31, 

5 

5.44, 

6.00, 

6.62, 

7.32, 

8.20, 

9.75, 

11.02, 

10.70, 

10.13, 

P.45, 

8.70 

* 

’  .72, 

7.11, 

6.40  , 

7.  4 S  , 

3.7<s, 

4.20  , 

4.  72, 

5.  50, 

6.70, 

7.82, 

7.  35, 

6.  70, 

5.85, 

5.12 

• 

4. 4C  , 

3 .74, 

3.18, 

1 

1.  ’6, 

1.41, 

1.56, 

1.  79, 

2.  14, 

2.72, 

3.54, 

2.R2, 

2.  31, 

1.67, 

1 .2  6 

• 

.  °2, 

0.65, 

'’.53, 

20.,a.,o.,o.,o.,o.,o.,c.,a.,j..o.,').,o.,o.) 

CATA  ((  ( S ( J , I )  ,1  =  1,14)  ,  J=  1 , 5  )  = 

71  4.  61,  14.  7 1,14. 89, 15. 37,  15. ?4,  15.  4  2,  15.59,  15.  77, 15. 93, 16. 10,  16.  26 
*16.  41,16. 58, lb. 73, 

S14. 58, 14. 77,14.  <36, 18.1 3, 1<;.U, 15.  49, IS. 68, IS.  95, 16. 02, 16. 20, 16. 36 
*16. 31, 16.67 ,16.81, 

314. 66, 1^. 83, 15. 04, 15. 2  2,  lc. -*1,15.59,15.78, 15. 96, 16. 13, 16.. 30, lo.47 
*16.61,16.77,16.90, 

1 14.  74,  14. «4,  15.  14,  15. 33,  15. 53,  15. 72, IS. 92,  16.  11, 16. 29, 16. 44, 16. 60 
*16. ”3, 16.8’, 16.98, 

214. 79, 14. 99, 15. 19, 15. 40, 15. 6n ,15. 53, 16. 06, lo. ’2, lb. 37, 16. 50  ,16. 6* 
*16.76, 16.89,16. 99) 

I F ( NC A  Sfc . FQ . 2 )  GO  TO  20 
1C  T=(x-32.)*S./9.*273. 

11  Js  3 , 5 

Ip ( x . GE. 9  (  J) )  GO  TO  12 

11  CONTINUE 

12  nn  13  1=3,14 
IF(T.LE.TKII) )  GO  TO  14 

13  CONTINUE 

14  L?=T 

L 1= I -2 

no  15  L=Li,L2 

CALL  TRIQUA (B  ( J-2)  ,3(J-t),B(J),DTK(J-2,L),OTK(J-l,L),nTMJ,L),Al, 

*  A  5 ,  A3) 

CALL  triuuA(B(J-2) , a ( J- 1 ) ,6( J),S(J-?,L) ,S (J-l ,  L >  ,  S ( J, L) ,81,B2,BX) 
OTL  (L ) =A 1*X*X»A2*X^  A  3 
STL  (L ) =31*X*X*92*X»93 

15  CONTINUE 

CALL  TRIUUA  (TK(I-2)  ,  TK(T-l)  ,  T<(I),r>TL  (I  -2)  ,  OTL  ( 1-1)  ,  9TL  (  I)  ,C1,C2, 

*  C  3 ) 

CALL  TOIOIJA  (TK(I-2)  ,  TK(  1-1 )  ,  Tw(  I)  ,  SIL  (1-2)  ,SIL  (1-1)  ,SIL  (I ) ,  01, 02, 

*  C’l 

C T  =  C.  1  *  T *  T  ♦  ,?*T4C3 

0T=0T*9./5. 

ST=01*T*T02*T*03 

RFTUPN 


20  T=(Y-32. )*5./9.*273. 

DO  21  1=3,14 

TF (T.LE.TK(I) )  GO  TO  22 

21  CONTINUE 


Copy  available  (o  DTiC  w.. 
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22  00  2 ?  J=l,5 

CALL  TRIQUA (TK(I-2)  , T K  C I -1 > , TK(I),DTK (J ♦  1-2)  ,  QTK(  J,  1-1) «OTK(J*I)  . 
#  El »  £2 »  E3 ) 

CALL  TRIQUA  (TK (1-2) ,TK(I-1) ,  TK  (I )  ,  S  (  J,  I -2) , S( J, 1-1) ,S (J, I) , Fl,F2, 
*F3) 

OTL ( J)=E1*T*T+E2»T*E3 
SIL  (J)=Fl*T*7>F2*T>E3 

23  CONTINUE 

00  24  J=3 » 5 

IF(X.LE.SIL (J) )  GO  TO  25 

24  CONTINUE 

2  5  CALL  TRIQUA (SIL (J-2) ♦ SIL ( J-i > , SI L( J) ,DTL(J-21  , OTL  ( J-l ) , OTL ( J) , 
*G1«G2«G3) 

CALL  TRIQUA (SIL ( J-2 ) »  SIL ( 2-1) »SIL( J)  tB( J-2)  ,8 (J-l) , 3 ( Jl  , H 1 , H 2 , H3 ) 

DT=G1*X*X+G2*X*G3 

0T=0T*9./5. 

3T=Hl*X*X+H2*y*H3 

o" TURN 

END 


SUBROUTINE  T RI QUA { X 1 , X2 , X3,Y1,Y2,Y3,A,B,C) 

THIS  SUBROUTINE  FINOS  COEFFICENTS  FOR  A  QUADRATIC  EQUATION 
A  =  ( 1./ (X2-X3) ) * ( (Y1-Y2) / ( X 1-X  2) - ( Y1-Y3 ) / ( XI- X3 )) 

B= (Y1-Y7) /(Xl-X3)-A*  (X 1 ♦ X3) 

C=  Y2-A*X2*X?-3*X2 

RETURN 

END 


Copy  available  to  D'l'iC  does  not 

permit  fully  legible  reproduction 
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SUBROUTINE  GRNBRO 

-  |  C  ™ I S  SUBROUTINE  PLOTS  THE  OUTPUT  OP  THE  MAIN  PROGRAM. 

I  COMMON  OCX  (50  ,  <50  »  ,  VElX  (BO,  SOI  ,NRX,NV  (50  )  ,  COPX  (50, 50  ) 

COMMON  WC,LC,TC,TP,TFH,TFC,OTM,OTW,RCU  j 

COMMON  N , NN  ! 

CIMcNSION  VELXXC50) ,QCXX(5C) ,COPXX(50> , IPAK  (120 J  !■ 

REAL  LC  i 

CALL  OOMOPS  ! 

CALL  BGNPL ( N) 

CALL  PAGE <11. 

CALL  PHYSOP <9.8, 1 .5) 

CALL  N03R0° 

CALL  nANGL: <9n.) 

CALL  TITLE!"  " , 0 , "G A OOL I N I U“  PLATE  CHANNEL  VELOCITY  -  r T/SECS" , 1 00 
l, ‘TOOLING  CAPACITY  -  ETIJ/Hk/Ll  OP  GO  <SOLIO  L  T  NES  )  J"  ,  10  0 , 6.  ,  8  .  ) 

CALL  HEIGHT!. 13) 

>  CALL  MESSAG<"MAGN£TO  T  Hrrt  MA  L  HEAT  PUMP  PERFORMANCE  AT  ,  1  00  ,  .  2 , 9  .  >  j 

CALL  REALNO  (3CU,  1 ,"A  3UT"  ,"A  RUT")  jj 

;  CALL  M £ S S A G  < ”  TESLAS",  lJn,*'AjuT",,*Ar*UT”)  M 

CALL  HEIGHT  (  .P  9)  jl 

j  CALL  MESSAG! "CHANNEL  HEIGHT,  TN.  =  «",  100, -.  7,  8.  R)  j  i 

1  CALL  REALNO  <WC, 2, "ABUT", "AP'JT")  |! 

■  j  CALL  M ESS AG l"  °L  A  T£  SPACING, IN.  = ?“ , 10 1 , "A  RUT" , " AB UT" )  M 


CALL  REALN0!TC,3,"A3UT",“APIJT") 

CALL  h£SSAG<“  TEMP, HOT  ENO,F  =3”,  10 0 , “ AOU T“ , " ABUT" ) 
CALL  REALNO  <  TFH,  0  ,  "  A  8UT  "  *  "A  afJT"l 

CALL  MESSAG!"  MIXING  T^hP, F  =  $", l Q  0 ,“AbUT~,“A0UT") 

CALL  REALNO (OTM, 1 , " A3 UT" , " AO UT" ) 

CALL  MESSAG ("CHANNtL  LENGTH, IN, =$", 100»-.T*8,3) 

CALL  RFALN0!LC*2," ABUT"," ACM*") 

CALL  MESSAG!"  PLATE  TH IC KNE3S, IN .=*“, 1 00 , "A BUT" ," ABUT* 
CALL  REALN0!TP,3,"A8UT”,"A?Ut") 

CALL  MESSAG!"  TEMP  COLD  E  4C  ,F  = J" , ID  A , " A  BUT" ," AB UT”) 
CALL  REALNO(TFC,0,"A8UT","A  I  UT" I 

CALL  MESSAG!"  H . E . OEL T A -T , F  =  «” , 100 , " A BUT" ,” ABUT” ) 


HPL=  0.  1 

CALL  HEIGHT  t  HPL I 

inUMM Y=L I N£ ST !IPAK, 120, HP ) 

CALL  LINES ( 13H  1-CYCLE/MIN" ,  IPA<,1) 
CALL  LINES ( 13h  2“CYCLE/MIN„ , 1° AK,2 ) 
CALL  LINES(13H  4-C YCL E/ M INI  ,  I  =AK , 3 > 
CALL  LINESI13H  8-CYCLE/ MIN* , TPAX, 4) 
CALL  LINES!  13H12-CYCL E/M  IN  -  ,  IDAY,5) 
CALL  LINES (13H18-CYCLE/ MIN* , I  PAX, 61 
CALL  LINES! 13H24 -C YCL E/MI N*  ,TPAK,7) 
CALL  LINES!13H30 -CYCLE/MI NS ,IDAK,8) 
CALL  LINES (13H40-CYCLE/MTN* ,IPAK,9) 
CALL  LINFS  <13H90- CYCLE/ MIN* , TPAK,10> 
CALL  LINES ( 1 3H60 -CYCL E/MIN? ,IPAK, 11) 
XHsXLEGNO !IPAK,NRX*l) 

YH=YLEGNO !IPAK,NRX+1 ) 

XWA=5. 3-XW 
YHA=T.a-YH 

CALL  3LNK1!-7.9,-7.7hYH,S.7-XW,5.9,1> 
CALL  FRAME 

CALL  RESET <~HEIGHT") 
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CALL  XTICKSt  9  > 

CALL  YTICKS(5) 

CALL  XLOGt .3,3. ,0., 75.0) 

CALL  OOT 
CALL  GRI0(1,1> 

CALL  RES-T(-OOT“) 

CALL  SPLINE 
CALL  NOCHr  K 
DO  10  NC=1,NRX 
NPX  =  NV(  NCJ 
00  8  1=1, NPX 
VELXX(I)=V-.X(NC,I) 

QCXX(I)=QCX(NC,I) 

8  CONTINUE 

CALL  CURV£(VELXX,QCXX,NPX,3) 

10  CONTINUE 

CALL  YTICKS(?) 

CALL  YGRAXS(0.  , 1., 12. ,8. .-COEFFICIENT  OF  PERFORMANCE  (DASHEO  LINES 

1) 3”, -100, 6., 0. ) 

CALL  OASH 
DO  20  NC=1,NRX 
NPX=NV (NO) 

CALL  MARK-R(NC-l) 

00  18  1  =  1, NPX 
VELXX(I)=VF_X  (NC.I) 

COPXX  (I )=OOPX  (NC.I  ) 

18  CONTINUE 

CALL  CURVE  (V-LXX,  COPXX,  NPX,  2) 

20  CONTINUE 

CALL  RESET (-3LNK1”) 

CALL  HEIGHT(HP_) 

CALL  LEGENO(I3AK,NRX,XWA,YHA) 

CALL  ENPPl(N) 

IF(N.EQ.NN)  GO  TO  25 
RETURN 

25  CALL  OONEPL 
RETURN 
END 


^  does  ** 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1  OTNSROC  REPORTS.  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT 

2  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY.  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  COOE  OF  THE  ORIGINATING  DEPARTMENT  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS 
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